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THE PSEUDO-EVOLUTION OF HALO MASS 
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ABSTRACT 

A dark matter halo is commonly defined as a spherical overdensity of matter with respect to a reference 
density, such as the critical density or the mean matter density of the Universe. We show that such definitions 
lead to a spurious evolution in the halo's mass even if its physical density profile remains constant over time. 
This pseudo-evolution in mass is caused by the evolution of the reference density with redshift, and has little 
connection with the actual physical accretion of mass. We compute the pseudo-evolution of halos identified 
in a large N-body simulation from z - 1 to 0, and show that it increases halo masses significantly across a 
wide range of halo masses and overdensities. Pseudo-evolution accounts for almost the entire mass evolution 
of halos with M^oap ~ 10 12 /2 _1 M , while for larger halos it still accounts for ~ 50% of their overall mass 
evolution. We estimate the magnitude of the pseudo-evolution assuming that halo density profiles remain static 
in physical coordinates, and show that this simple model predicts the pseudo-evolution of simulated halos to a 
few percent accuracy. We discuss the impact of pseudo-evolution on the evolution of the halo mass function. 
We show that the non-evolution of the low-mass end of the halo mass function is the result of a fortuitous 
cancellation between pseudo-evolution and the absorption of small halos into larger hosts. We also show that 
the evolution of the low-mass end of the concentration-mass relation observed in simulations is almost entirely 
due to the pseudo-evolution of mass. Finally, we discuss the implications of our results for the interpretation 
of the evolution of various scaling relations between the observable properties of galaxies and galaxy clusters, 
and their halo masses. 

Subject headings: cosmology: theory - methods: numerical - dark matter - galaxies: halos 



1. INTRODUCTION 

In a cold dark matter cosm ological scenario (see e.g. 
lPeebleslfl98llDavis et al.ll 19851) . the drama of galaxy forma- 
tion unfolds at the virialized peaks of the density field, or ha- 
los. Although galaxies themselves are highly diverse, a num- 
ber of their properties exhibit remarkable regularity and can 
be expressed as galaxy scaling relations. In particular, the 
stellar mass-halo mass relation and the luminosity-halo mass 
relation of central galaxies constrain important aspects of 
galaxy formation and have be en studied via a va r iety of probes 
such as satellite kinema t ics dPrada et al.l 120031 : IConroy et alJ 
120071: lMore^taLl|2009l 1201 lbl) galaxy- galaxy weak lens- 
ing dSeliakl |200a iMcKav et al.l 1200 it Mandelbaum et alJ 
I2006t iParker et al.l 120071: ISchulz etal.1 l2010t). the abun- 
dance of ga laxies and their clustering (lYang et alJ 1 2003 1: 
Zehavi et all 12004 ITinker et al.l l200l IZehavi etaLl 120051: 
Skibba et all boOfi Ivan den Bosch et alJ 120071: iBrown et all 
2008t IConrov & Wechslerl 120091: iMoster et al.l 120101 1201 " 



such as X-ray temperature, entropy profile, the mass of the 
intracluster-gas, or their evolution w ith redshift, are often 
described using a self-sim ilar model dKaiserj[l986l see also 
iKravtsov & Borganl 120121 for a review). This model pro- 
vides predictions for the scaling relations between halo mass 
and the observable properties of clusters. Lar ge observational 
campaigns have been undertaken in the past (|Vikhlinin et alJ 



I2006t iBohringer et al J 120071: iMantz et al 



Benson et al.l 



201C ) and are also 



2011) to calibrate 



Beh roozi et al. 120101: | Yang etail l2012al) or a combination 
of the above probes ( Yooetal. I2006T Cacciato etal.1 120091 : 
iLeauthaud et al.ll2012t More et al.ll2012al) . In order to under- 



stand the formation and evolution of galaxies, it is crucial to 
interpret the evolution of these scaling relations, which in turn 
requires a solid understanding of the evolution of halo masses 
with cosmic time. 

Analogously, the largest halos in the Universe host clus- 
ters of galaxies, which themselves serve as laboratories for 
galaxy formation. The observable properties of clusters, 
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currently under way (e.g. 
these scaling relations since they are necessary to obtain cos- 
mological constraints from the ob served abundance of c lus- 
ters and its redsh ift evolution (e.g.. lVikhlinin et aTll2009aL see 
I Allen et al.l201 ll for a recent review). However, such observa- 
tional campaigns must be supplemented by sound theoretical 
models for the evolution of the scaling relations, which have 
still not been f ully developed (see, e.g., a recent analysis by 
iLin et al.ll2012h . 

When quoting the scaling relations between halo mass and 
galaxy (or galaxy cluster) properties, observers inevitably 
adopt a specific definition for the boundaries of halos, often 
based on the extent of their observations. However, numerical 
simulations show that dark matter halos exhibit smooth den- 
sity profiles without well-defined boundaries, which makes 
the definition of the halo boundary and the associated halo 
mass ambiguous. The mass definition often used in the litera- 
ture corresponds to the mass within a spherical boundary that 
encloses a given overdens ity, A(z), with respec t to a reference 
density, p re f(z) (see, e.g., ICole & Lacey||1996l) . This spheri- 
cal overdensity (SO) halo mass, Ma(z), and radius, Ra(z), are 
thus related via the following equation, 



M A (z) 



jTRl(z)A(z)PrAz). 



(1) 



The most common choices of reference density are either the 
critical density, p c , or the mean matter density, p, of the uni- 
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verse at a given cosmic epoch. The parameter A can be chosen 
arbitrarily, but certain values such as A = 180 can be justi- 
fied with the spheric al top hat collapse m odel for an Einstein- 
de Sitter cosmology dGunn & Gottlll972h . The spherical col- 
lapse model has also been generalized to cosmological models 
which include a cosmological constant or non-zero curv ature 
(lLahavet all 19911: lLacev & Cold[l993l: lEke et al.l[l996l) . 

The fundamental issue with the mass definition of Equa- 
tion ((TJ is that the reference density evolves with cosmic time, 
leading to an evolution in halo mass even if the physical den- 
sity profile of the halo is constant. For the remainder of this 
paper, we shall call the evolution of halo mass due to changing 
reference density pseudo-evolution because it is due solely to 
the mass definition and not to any actual physical mass evo- 
lution caused by the accretion of new material. Note that the 
actual evolution of spherical overdensity mass, which we shall 
call mass evolution, is a combination of the physical evolution 
due to the accretion of matter and pseudo-evolution (see Ap- 
pendix [A}. 

The fact that the evolution of the SO mass may not cor- 
respond to any actu al physical evo l ution of mass has been 
pointed out before. iDiemand et al.l (120071) analyzed the ac- 
cretion history of the Milky Way sized Via Lactea halo and 
found no significant physical growth after z = 1, even 
though the virial mass of the halo increased significantly. 
iPrada et ail ( 120061) studied the outer regions of collapsed ha- 
los at z = and found no systematic infall for halos with 
masses lower than 5 x 10 12 h~ l M Q . In a follow-up study, 
ICuesta et al.l (120081) demonstrated a lack of physical accretion 
onto galaxy mass halos, and proposed an alternative mass def- 
inition that aims to include all m ass bound to a halo (see also 
lAnderh alden & Diemand 120111) . Although such a mass def- 
inition may be more physical and closer to the meaning of 
mass in analytical models of halo collapse and evolution, its 
observational analogue is very difficult or even impossible to 
measure for real systems. Thus, the SO mass is most often 
used in observations, and a proper interpretation of observa- 
tional results should take into account the pseudo-evolution 
inherent i n this mass definition . For th e case of cluster scaling 
relations. iRravtsov & Borganil(l2012|) argued that part of their 
evolution is due to pseudo-evolution. 

In this paper we seek to quantify the pseudo-evolution 
of the SO halo mass due to an evolving reference den- 
sity. In particular, we focus on quantifying the contribu- 
tion of the pseudo-evolution to the mass accretion history 
(MAH) of halos, since the MAH describes the dark matter 
part of the observable mass scaling relations for galaxies and 
galaxy clusters. Much work has been invested into quanti- 
fying MAHs, but the contributions from physical accretion 
and pseudo-evolution are generally not dist i nguished (see e.g . 
Wechsler et aT]|2002T: Ivan den Boschl l2002T: IZhao et all 12001 
Miller et al.ll2006t IZhao et al.ll2009ir We will also investigate 
the impact of pseudo-evolution on the evolution of the halo 
mass function. 

This paper is structured as follows. In §|2] we describe a 
simple analytical mo del based on the NFW density profiles 
dNavarro et al.l IT997) to gauge the contribution of pseudo- 
evolution to the total MAH. In §0 we quantify the pseudo- 
evolution of actual halos in cosmological simulations and 
show that the results are in good agreement with the analytic 
model of §|2] We discuss caveats and implications of our re- 
sults as well as directions for future work in $4] an d gi ye a 
summary of our results in §0 Throughout this paper, we de- 
note overdensities as A c if they are defined relative to p c , and 
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Fig. 1 . — Visualization of the static halo model. The solid line shows the 
spherical mass profile as a function of enclosed density (M/V) for a halo of 
mass M200p = 2 X 10 12 /!~'M Q from the Bolshoi simulation at z = 0. The x- 
axis is reversed, so that the left side of the plot corresponds to the high-density 
center of the halo, and the right to the low-density outskirts. The vertical 
dashed lines indicate the spherical overdensity 200p at redshifts 0, 0.5 and 
1, and the horizontal dashed lines mark the corresponding halo mass M200p- 
Since the reference density p oc ( 1 + z) 3 , the halo density threshold increases 
with redshift, and Mioojs decreases. Even if the physical mass distribution 
of this halo was kept fixed between z = 1 and z = 0, its mass M200P would 
undergo a pseudo-evolution from 9.6 X 10 1 1 /r'M to 2 X 10 12 fc _1 M s . 

A m if they are defined relative to p. We also use A v i r to de- 
note the redshift and cosmology-dependent virial overdensity 
predicted by the spherical collapse model, which corresponds 
to A v j r * 358 at z — and A v ; r =s 180 at z > 2 with respect 
to the mean background density for the concordance fiducia l 
cosmology used in this paper (e.g. jBrvan & Norman 19981) . 
All densities and radii are expressed in physical units, unless 
stated otherwise. 

2. THE STATIC HALO MODEL 

Let us consider a density peak in the Universe around which 
the matter density profile in physical units has not evolved 
since a given redshift, zo- As shown in Fig.[T| the halo mass 
associated with this density peak will change purely due to the 
evolution of the reference density used to define its boundary. 
This evolution in mass can be quantified using the density pro- 
file of the halo at redshift zo- Let us assume that the density 
distribution around this density peak is described by the uni- 
versal density profile given by 

p(r ' Zo) = < / uf / ^ ' (2) 
(r/r s )(l +r/r s Y 

which has been found to be a reasonable approximation of the 
typical density pr ofiles around density peaks in cold dark mat- 
ter cosmologies dNavarro etal J [19971 NFW hereafter). The 
scale radius r s and the halo radius are related by the con- 
centration parameter ca = Rhlr?,- Under our assumption that 
the density profile around this peak does not evolve and that 
profile of Equation d2} is a good description of the actual pro- 
file at the radii of interest, the halo mass at any redshift z can 
be expressed in terms of the characteristic density p s and r s , 
by integrating this static density profile within the halo radius 
Ra(z), such that 

M A (z)= p(r)4jrr 2 dr = p s 4^p[c A (z)] , (3) 

Jo 
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Fig. 2. — Predictions of the static halo model. Lines show the pseudo- 
evolution of halo mass due to changing reference density as a function of 
E(z) relative to the halo mass at z = 0. The labels indicate the halo mass in 
log( /r'M Q ). From top to bottom, the panels show the evolution of Mioop, 
M5oo Pc and A/ V j r , where p and c indicate the average and critical densities of 
the universe, respectively, and M v j r c orresponds to an evolving overdensity 
according to Bryan & Norman 1 1998 ) . For z < 0.5, the pseudo-evolution is 
largest for Mioop, but the overall trend is the same for all three mass defini- 
tions. 

where the function //[x] is given by 

/iW = ]n(l+jr)--^-. (4) 
1 + x 

Equating the right hand sides of Equations (Q]) and (0, we 
obtain a relation between the concentration parameter of the 
halo at redshift z and the concentration parameter at redshift 

Zo, 



3p. 5 



ca(z) 3 



//[c A (z)] Mz)f>mf(z) 

ca(zo) 3 



M[ca(zo)] 



A(Zo)Pref(Zo) 



Mz)p ref (z) 



(5) 
(6) 



This relation can in turn be used to find the evolution of halo 
mass according to the equation 



M A (z) = M A (zo) 



/4£a(z)] 
m[ca(zo)] 



(7) 



As examples, we consider three commonly used definitions 
of halo mass in the literature, (i) A m (z) = 200 as in studies of 
the halo occupation distribution of galaxies, (ii) A c (z) = 500 
as in studies involving galaxy cluster observations, and (iii) 
A c (z) = Avir . We assume that the density profile around a den- 
sity peak does not change from redshift z = 1 to z = 0. With- 



out loss of generality, we use zo — to define the static density 
profile in physical units. We consider t he con centration-mass 
relation at z — given by IZhao et a D d2009l hereafter Z09), 
and use Equations (O and © to obtain the mass evolution 
due to pseudo-evolution. 

Fig. |2] shows the mass evolution histories for halos of dif- 
ferent masses as predicted by our static halo model. The dif- 
ferent panels correspond to the three commonly used over- 
density definitions. Each panel shows the pseudo-evolution 
of mass from z = to z — 5, normalized to the halo's mass 
at z = 0, as a function of the expansion rate in units of the 
Hubble constant, E{z). Assuming a flat ACDM cosmology, 
E(z) is defined as usual, 



e(z) = Vn A + fi„,(i +z) 3 



(8) 



As expected, more massive halos undergo a larger evolution 
due to the lower values of their concentrations. Regardless of 
the exact mass definition, the fractional change in halo mass 
due to pseudo-evolution can be as large as ~ 0.5 by z = 1 . 
The shape of the mass evolution history with redshift is not 
only a function of halo mass, but also depends upon the exact 
mass definition. Its functional form is better approximated by 
a power law of E{z) than (1 + z), but still shows deviations 
from an exact power law behaviour. 

3. HALO MASS EVOLUTION IN SIMULATIONS 
3.1. Numerical simulation 

To quantify the pseudo-evolution of mass using realistic 
halo profiles, we use a sample of halos extracted from a dis- 
sipationless cosmological simulation of the A CDM model. 
Speci fically, we use the Bolshoi simulation (IKlvpin et al.l 
1201 ll) . which followed the evolution of the matter distri- 
bution using th e Ad aptive Refinement Tree (A RT) code 
dKravtsov et al.l 119971: iGottloeber & Klvpinl 120081) in a flat 
ACDM model with parameters Q m = 1 - Qa = 
0.27, £l b = 0.0469, h = H /(100 kms" 1 Mpc" 1 ) = 
0.7, erg = 0.82 and « s = 0.95. These cosmolog- 
ical para meters are compati ble with measurements from 
WMAP7 (iJarosik et al.l 120111) a com bination of WMAP5, 
BA O and SNe (iDunklev et all 120091) X-Rav cluster stud- 
ies dMantz et al.ll2008t IVikhlinin et al ll2009bl) . and observa- 
tions of the clustering of galaxies and galaxy-galaxy/cluste r 
weak lensing (see e.g. lTinker et alJ|2012tlMore et alj|2012bl) . 
The same cosmology was used for the calculations shown in 
Fig. 12 and will be used for the remainder of this paper. 

The Bolshoi simulation uses 2048 3 ~ 8 billion particles 
to follow the evolution of the matter distribution in a cubic 
box of size 250 /i~'Mpc on each side, which corresponds to 
a particle mass of 1.35 x 10 s /z _1 M . This implies that the 
smallest halos considered in this paper (M v ; r = 2xlO 11 /i~ 1 M ) 
are resolved by over 1000 particles. As density peaks, we use 
the centers of halos from a catalo g generated with the bound 
density maxima (B DM) algorithm (IKlvpin & Holtzmanl 19971 : 
IKlvpin et al.ll201lb . 

We identified all distinct halos with M vil > 2x 10 1 1 M Q from 
the simulation at z = 0, resulting in a sample of about 240, 000 
halos. For each of these halos, we constructed radial density 
profiles by summing the particle contributions in 80 logarith- 
mically spaced bins, spanning radii from 0.05 to 10 R v i T . As a 
sanity check, we extracted density profile s for the same Bol - 
shoi halos from the MultiDARK database dRiebe et alj|2011l) . 
compared them to ours and found excellent agreement. The 
larger radial range and finer resolution of our profiles com- 
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pared to the ones existing on the database allowed us to define 
halo masses for lower overdensity thresholds of A m = 200, 
and much lower background densities than those used to de- 
fine the catalog. 

For the purposes of this paper, we focus on the mass evolu- 
tion from z = 1 to z = 0. Before z ~ 1, £2 m ~ 1, meaning that 
physical accretion dominates for halos in the mass range of in- 
terest (Z09). Thus, the impact of pseudo-evolution is expected 
to be largest from z = 1 to z = 0. Furthermore, most obser- 
vational determinations of galaxy and galaxy cluster scaling 
relations have been performed below z = 1, which makes this 
particular redshift range the most interesting to investigate. 

3.2. The evolution of M200P and M^qo Pc 

In a first experiment, we focus on mass definitions based 
on the mean matter density of the universe, such as M2mp- 
To estimate the pseudo-evolution of halo mass using realistic 
matter density profiles, we assume that the density profiles 
stay constant in physical units from z = to z — 1 and evolve 
the background density according to 

p(z) = (l+z) 3 p(0). (9) 

The radius, /?200p, of the halo is then numerically identified to 
be the radius which encloses an average overdensity of A = 
200 with respect to p(z), and the mass, M2oop, follows from 
Equation (fTJ. 

In the left hand panel of Fig. [3] we show the mean of the 
ratio of the evolved halo mass to the halo mass at redshift 
zero as a function of redshift using open circles, and the 16 
and 84 percentiles of the distribution of this ratio as the er- 
ror bars. The analytical estimate from §|2]is shown using a 
solid line, with gray contours indicating the 68% confidence 
interval. This uncertainty in the analytical prediction is due to 
scatter in the concentration-mass relat ion, which we assume d 
to be 0.14 dex based on the results of IWechsler et al.1 (120021) . 
The excellent agreement between the analytical estimate and 
the results from the halo profiles implies that our assumption 
about the density distribution at redshift zero (from the models 
of Z09) is not too far off from the actual density distribution of 
halos in the simulation. This shows that the analytical model 
can provide an excellent description of the mass evolution if 
the physical density distribution of the halos was indeed con- 
stant. 

Next, we would like to contrast the predictions of the static 
halo model with the true mass evolution histories of halos 
observed in simulations. In principle, we could derive the 
true mass evolution history of a halo by constructing merger 
trees from the simulation, and following the main progeni- 
tor branch of the halo. Instead, for simplicity, we make use 
of the Z09 model for the mass evolution histories of halos. 
This model has been shown to accurately reproduce the mass 
evolution histories for a large variety of cosmological models 
(scale-free or ACDM). The results from this model will thus 
include the effects of both pseudo-evolution and the actual 
physical accretion of mass. By comparing the mass evolution 
histories from our static halo model to these realistic mass 
evolution histories, we can disentangle the two effects. 

The mass evolution histories predicted by the model of Z09 
are shown by dashed lines in Fig. [3] The comparison clearly 
shows that a significant fraction of the halo mass evolution is 
due to the evolving reference density in the mass definition. 
In fact, for halos of mass M200/3 i 10 12 /i _1 M Q , the mass evo- 
lution history at z < 1 can be explained almost entirely by 



the static halo model, and true physical accretion constitutes 
a negligible portion of the overall evolution. For high-mass 
halos (M200P ~ 10 14 /i^'Mq), the pseudo-evolution accounts 
for roughly 60 percent of the mass evolution. 

There is, however, one potential issue with these results. By 
defining the static density profiles at redshift zero and work- 
ing backwards in time, we always shift the boundary of ha- 
los radially inwards to define halos at higher redshifts. This 
means that our experiment does not exclude the possibility 
that the density profile in the outer parts of the halo (beyond 
the boundary) change at higher redshifts. This could very 
well be the case, if halos underwent an inside-out growth by 
primarily accreting mass at the outskirts as time progresses. 
However, this issue can be easily addressed by defining the 
static density profiles at redshift zo — 1 and working forward 
in time. This allows us to test whether the density profiles of 
the halos follow the NFW form given by Equation (fJJ even in 
the outskirts of halos. 

The open circles in the right hand panels of Fig. [3] show 
the result of using the density profiles from the simulation at 
redshift zo — 1 and determining how the mass should grow 
as redshift decreases. We call this experiment a test of the 
forward evolution of the halo masses. Note that the forward 
evolution panels show the inverse of the evolved mass in units 
of the mass at zo, such that equal numbers in the left and right 
hand panels correspond to equal amounts of halo growth. As 
in the left hand panel, we show the prediction from our static 
halo model as a solid line. If the physical density profile in 
the outskirts of halos at z = 1 was different from what it is 
at z = 0, we would see a difference between the prediction 
of our analytical model (which assumes that the profiles are 
NFW all throughout) and the results from the actual density 
profiles at zo = 1 from the simulation. However, the excel- 
lent agreement between the two results shows that the density 
profiles indeed are described, to a good approximation, by the 
NFW form even in the outskirts of halos that we considered 
at Zo = 1. In Fig. [3] we also show the prediction of the mass 
evolution history from the models of Z09 using dashed lines. 
The fair agreement between our predictions of the static halo 
model with the mass evolution histories from Z09 at low halo 
masses (M200P ~ 10 12 h~ l M G ) validate our basic premise that 
the density profiles of such halos are already established at 
Z = 1 and do not change with time. 

The mass evolution shown in Fig. prefers to a sample of 
isolated halos only. Inevitably, this sample contains close 
pairs of halos which are not identified as overlapping in the 
halo catalog, but whose density profiles pick up contributions 
from particles belonging to the other halo. As the density 
profiles we use extend as far as 10 R w i r , this is the case for a 
significant number of halos in our sample. When we evolve 
such a halo forward in time, its radius grows, and the con- 
tribution from other halos leads to excess mass growth, which 
manifests itself as a large scatter in the mass evolution history. 
However, many of these halos end up lying inside the radius 
of a neighbouring, larger halo, and we need to exclude them 
from the averaged mass evolutions shown in Fig. [3] 

In principle, the most straightforward way to identify sub- 
halos would be to use merger trees of the Bolshoi simula- 
tion. However, merger trees are generated with knowledge of 
the full mass evolution of halos (rather than just the pseudo- 
evolution), as well as their motion and acceleration. In the 
spirit of our extremely simple model of static halo profiles, 
we wish to avoid using such information, and rely only on the 
density profiles and initial positions of the halos in our sam- 
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Fig. 3. — Mass evolution for different halo masses for the M200/5 mass definition, in five logarithmic mass bins. The dashed lines show the actual mass evolution 
of halos as predicted by the Z09 model. The solid lines show the evolution of mass due solely to the evolution of the reference density in the mass definition 
(pseudo-evolution) as predicted by the static halo evolution model (see ^2). The gray band around the solid lines shows the 68% confidence interval due to 
scatter in the concentration-mass relation. The red points show the pseudo-evolution computed using density profiles from the Bolshoi simulation, with error bars 
indicating the mass range containing 68% of the halos in a given mass bin. The left panels show this evolution computed by extrapolating the mass evolution 
using profiles at zo = and going backwards in time, while the right panels show the evolution of profiles at zo = 1 . going forward in time. Note that the right 
hand panels show the inverse of M(z), meaning that equal values at z = 1 in the left hand panels and z = in the right hand panels indicate the same amount of 
evolution between z = 1 and 0. However, the mass bins on the left and right do not correspond to the same halo masses, and should thus not be compared directly. 
These results demonstrate that the pseudo-evolution accounts for about a factor of two at all hal o ma sses. For small halos (MiQQp i. 10 1 " ft Mq), this represents 
the majority of their overall evolution, while for very large halos it accounts for about half. See 33.2l for further discussion. 
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Fig. 4. — Like Fig. [3] but for the M5oo Pl mass definition. Only the highest and lowest of five mass bins are shown. The analytical prediction of the static halo 
model and the results from simulated halos match even better than for Mioop, and the Bolshoi results exhibit smaller scatter. This is caused by a smaller virial 
radius Rsqo Pc , meaning that irregularities in the outskirts of halos play a lesser role. The pseudo-evolution is a little weaker in Msoo Pc than in M200p, but still 
accounts for almost all of the mass evolution at low halo masses. 
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pie. 

In the case of backward evolution, halo radii decrease, 
meaning that we do not need to worry about halos becoming 
subhalos. In the case of forward evolution, we identify sub- 
halos as follows. At each redshift, we evolve the virial radii 
of all halos to match the evolved reference density. We then 
check whether the new virial radius of such a halo encloses the 
center of any other, smaller halos. Note that we assume that 
the neighboring halos stay at a constant physical distance and 
are not part of the Hubble flow, consistent with our assump- 
tion that the physical density around the peak does not evolve. 
We exclude the subhalos discovered in this manner from the 
current redshift bin, and all subsequent smaller redshifts. We 
start this process with the largest halo in the sample, marching 
down to the smallest halos. Once a halo has been found to be 
a subhalo, it cannot itself be the host of another halo. By the 
end of the evolution from zo = 1 to z = 0, a total of about 14% 
of the halos in the sample had become subhalos and been re- 
moved from the sample. The mass evolution averages shown 
in Fig. [3] are insensitive to the removal of subhalos, but the 
scatter is reduced significantly by this procedure. 

Following our discussion of mass definitions based on the 
mean matter density of the universe such as Mooop, we now 
investigate definitions based on the critical density. Note that 
the difference between those definitions is not only due to the 
different values for A which are typically chosen, but that p c 
evolves qualitatively different from p such that 



Pc(z) = P(Z) 



E\z) 



(10) 



Fig. [4] shows the evolution histories of Msoq Pc for two mass 
bins, compared to the true evolution represented by the model 
presented by Z09. The results were derived in exactly the 
same way as the results for M2oop, except for the different 
evolution of Ap re f. The pseudo-evolution is slightly weaker 
in Msoo Pc than Majop- This can be seen by comparing, for ex- 
ample, the lowest mass bins in Figs. [3] and [4] The weaker 
evolution in Msoop c may seem slightly counter-intuitive at 
first, since /?5oo Pl is smaller than 7?2oop (which implies that 
C500p t < QoopX an d mass profiles as a function of enclosed 
density tend to steepen towards the center of halos (see Fig.[T|i. 
However, the weaker evolution of p c compared to p more than 
offsets this effect. For example, at z — 1 the mean matter den- 
sity of the universe was a factor of 8 higher than today, but the 
critical density was only larger by a factor of E 2 (l) =s 2.9. 

3.3. The mass (non-) evolution of low-mass halos 

One of the most striking consequences of the results pre- 
sented in Figs. [3] and |4] is that the physical density profiles 
of low-mass halos (M200P ~ 10 12 /i~'M ) barely change after 
2=1, implying that they undergo negligible physical accre- 
tion. We seek to demonstrate this directly in Fig. [5] which 
shows the average radial velocities of particles around their 
halo centers at z — 1 . The colored lines correspond to three 
narrow mass bins. We extracted the velocity information from 
the Bolshoi simulation using the sam e binni ng scheme as for 
the density profiles. As ICuesta et al.l (120081) pointed out, the 
average radial velocities of low-mass halos amount to only 
a small fraction of V200 (the halo circular velocity at /?200p), 
while they grow to about half of V200 for the most massive bin 
in Fig.0 

We estimated the average mass infall rates at every radius 
using the density profile and the average infall velocity profile, 
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Fig. 5. — Mass infall around halos at z = 1. Top panel: Average radial 
velocity profiles of halos in three narrow mass bins, in units of the circular 
velocity V200 = VGM200/S200 (compare to Fig. 13 in Cuesta et al. 2008). For 
low-mass halos, the average infall velocity never exceeds a small fraction of 
their circular velocity. Beyond a few virial radii, the Hubble flow dominates 
and drags matter away from the halo. Center panel: average mass flux as 
derived from the velocities in the top panel, and density profiles. The out- 
flowing parts of the profiles (positive radial velocity) are omitted. Bottom 
panel: an estimate of the mass accreted since z = 1, We assume that the mass 
flux has remained constant since z = 1 (or 7.84 Gyr), and compute the total 
mass passing through a radial shell as a fraction of the original halo mass. 
The two arrows in the botto m panel indicate i?200p at z = (right arrow) and 
z = 1 (left arrow). See j|3.3l for a discussion of these results. 



v r (r), using 



dM 
dt 



Anr 2 p(r)v T {r) . 



(11) 



These mass infall rates are shown in the center panel of Fig. [5] 
for halos corresponding to the three mass bins. We estimate 
the total accreted mass since z — 1 by assuming that the ra- 
dial velocity profiles stay constant, the same way the density 
profiles remain fixed. Given this fixed inflow profile, we ex- 
trapolate the total mass passed through a radial shell by mul- 
tiplying the mass infall rate by the time since 2=1 (bottom 
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panel of Fig. [5j. 

We can now compare this estimate of the physically ac- 
creted mass to the mass evolution shown in Fig. [3] Accord- 
ing to Appendix [A] we need to consider the physical mass 
added within the radius the halo reaches at the end of its 
pseudo-evolution (see the second integral in Equation IA2b . 
For the case of backward evolution, this radius is Rwopiz =1) 
(marked with the left arrow in the bottom panel of Fig. 0, 
and we consider the infall rates at that radius. Note that the 
accreted mass is shown in units of the halo mass at z — 1. 
Since halos in the lowest mass bin, M200P ~ 2 x 10 11 /i _1 M , 
grow by about a factor of 2 by z — (mostly through pseudo- 
evolution), the infall rate at Rwop(z — 1) corresponds about 
10% of the halo mass today. This is in good agreement 
with the mass evolution results in Fig. [3] In the case of for- 
ward evolution, we need to consider the mass accreted within 
Rioopiz = 0) (indicated by the right arrow in the bottom panel 
of Fig. [5]). It appears that the accreted mass at that radius 
slightly overestimates the actual mass evolution for low-mass 
halos, and slightly underestimates the mass evolution of high- 
mass halos, which is presumably due to our assumption of 
fixed mass infall profiles. We intend to investigate this as- 
sumption in future work. 

3.4. Halo mass function from the static halo model 

Given the success of the static halo model in reproducing 
mass evolution histories for small mass halos, we would like 
to investigate the impact of the mass definition on the evo- 
lution of the halo mass function. For this purp ose, we start 
from t he zo — halo mass function calibration of iTinker et al.l 
J20081 hereafter T08), and use the mass evolution history 
(backward evolution) inferred from our static halo model to 
predict the resulting evolution of the halo mass function. In 
the static halo model, the mass assigned to density peaks be- 
comes smaller as redshift increases, which results in a mass- 
dependent shift of the halo mass function (MdN/dM) towards 
the left. In order to quantify this effect, we used the same pro- 
cedure as for the halo mass evolution. We selected the density 
profiles of halos of mass M v ; r > 2xlO 11 /?~ 1 M extracted from 
the simulation and calculated the expected evolution of mass 
assuming that the density profiles around peaks stay constant 
in physical units. 

Before comparing these results to the actual physical evolu- 
tion of the mass function observed in numerical simulations, 
we first establish that these results match the analytical pre- 
diction of Equation @. Since the corresponding mean mass 
evolution histories agree to a few percent (Fig. [3), we naively 
expect good agreement between the mass functions as well. 
However, Fig.[3]also reveals significant scatter in the pseudo- 
evolution of simulated halos, which could cause disagreement 
with the analytically predicted mass function. The compari- 
son between the two estimates is shown in Fig. [6] The an- 
alytical prediction of the static halo model for z - 0.5 and 
Z = 1 is shown with solid lines, while the actual halo mass 
function inferred from the pseudo-evolution of profiles in the 
simulation is shown by points. The fractional residuals be- 
tween the evolution from the actual density profiles and the 
prediction from the static halo model are shown in the bottom 
panel. At both redshifts, the prediction from our static halo 
model agrees extremely well (to better than 5 - 10%) with 
the expected evolution from the actual profiles of the halos at 
z — 0. For the case of mass definitions using higher over- 
densities, such as Msm Pc , we expect even better agreement, 
since the static halo describes the pseudo-evolution of halos 



in simulations more accurately (see Fig.|4]i. 

The right hand panel of Fig. [6]shows the same comparison 
as the left hand panel, but for the case of forward evolution. 
For the prediction of the static halo model, we start from the 
halo mass function of T08 at zo = 1, and evolve forward in 
time to z — 0.5 and z — 0. The slightly larger discrepancy at 
low masses at z = (~ 10%) is an artefact due to the mass 
cut of M > 2 x 10 11 /t'Mq that we imposed at zo = 1. Some 
halos which were excluded due to their low mass would have 
developed masses included in the shown range, leaving the 
low-mass end depleted. 

Note that, unrelated to the static halo model, there is a dis- 
crepancy of ~ 5 - 10% between the zo = and zo - I calibra- 
tion of the T08 mass function and the mass function obtained 
from the Bolshoi simulation. Given this initial offset at zo, we 
cannot expect a smaller discrepancy at subsequent redshifts. 

Furthermore, we considered the impact of statistical bias 
due to the presen ce of scatter in the mass evolution histories 
dEddingtonl 1 1 9 1 3b . It is evident from Fig. [3] that there is a 
non-negligible scatter in the mass evolution histories, and that 
this scatter is somewhat larger in the forward evolution case 
than in the backward evolution case. As the number density 
of halos is a decreasing function of halo mass, the number of 
halos that get up-scattered into a particular mass bin can be 
larger than the number of halos that get down-scattered out 
of that mass bin. However, we found that this bias does not 
influence the results appreciably. 

3.5. Comparison with the true mass function 

Having convinced ourselves that the halo mass functions 
predicted by the static halo model and halo profiles are con- 
sistent, we now wish to investigate the impact of pseudo- 
evolution on the true evolution of the halo mass function, 
analogous to the comparison of our mass evolution history 
with that predicted by the models of Z09. We use the cali- 
bration of the mass function provided by T08 to reflect the 
true evolution of the mass function measured in simulations. 
In Fig. [7] we present the comparison of the T08 mass func- 
tion at three different redshifts with the evolution of the mass 
function due to the pseudo-evolution of density profiles from 
the simulation. Note that the residual panel includes a wider 
range than in Fig. [6] Let us first focus on the left hand 
panel which shows the backward evolution case. Since our 
estimates of the mass evolution histories of low mass halos 
matched those observed in simulations (Fig. |3), we expect 
good agreement with the mass function at the low-mass end, 
and discrepancies predominantly at the high-mass end. The 
left hand panel of Fig. [7] however, reveals significant dis- 
agreement at both mass ends. Since the cause of deviations 
are different at the low and high-mass ends, we shall discuss 
those regions separately. 

At the high-mass end, we have shown that the growth of 
high-mass halos is largely due to physical accretion rather 
than pseudo-evolution. This rapid growth implies that the 
T08 mass function decreases strongly with redshift. As we 
evolve backwards in time, the progenitors of high-mass halos 
would need to be more massive than in reality if the growth 
was solely due to changing mass definition. The static halo 
model, therefore, overpredicts the number of large halos at 
z= 1. 

At the low-mass end, computing only the pseudo-evolution 
underpredicts the T08 calibration at z = 0.5 and 1 by roughly 
20 - 30% (left panel of Fig. |7J. From the calibration of 
T08, it appears that the number density of low mass halos 
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(M < 2 x 10 12 /i~'M ) stays constant since z = 1. The com- 
mon explanation for this observed non-evolution is that low- 
mass halos have already collapsed and do not physically grow 
in mass or number. However, it is clear from Fig.[3]that these 
halos do indeed undergo a significant mass evolution just due 
to pseudo-evolution. This conflict can be resolved by noting 
that, in the backward evolution case, the radii of the halos re- 
duce as we evolve the masses to higher redshifts. This may 
uncover substructures in the outer part of the halos, which 
can potentially be counted as isolated halos at those higher 
redshifts. Thus, the non-evolution of the halo mass function 
at the low mass end must be a result of the fortuitous cancel- 
lation of the effect of mass evolution and the addition of low 
mass halos at the outskirts of bigger halos. Since we limited 
our analysis to use only isolated halos at zo = 0, we were not 
able to quantify this effect in the case of backward evolution. 

However, we can investigate this effect by taking the static 
density profiles at zo = 1 and evolving them forward to z = 0.5 
and (forward evolution). In this case, the low-mass halos in 
the outskirts of larger mass halos should get absorbed. When 
we calculate the masses of halos by using static density pro- 
files, we partially account for this effect by removing small 
mass halos whose centers end up within the radius of larger 
mass halos at z = 0.5 and as discussed in j j3.2| The right 
hand panel of Fig. [7] shows the comparison between the for- 
ward evolution of the mass function of halos from the Bol- 
shoi simulation, and the T08 mass function. The compari- 
son shows that the discrepancy at the low mass end still per- 
sists, even after removing low-mass subhalos. The reason for 
this discrepancy is our implicit assumption that these density 
peaks are stationary. While we remove some subhalos as they 
are absorbed into larger halos, more halos which appear iso- 
lated at z — 1 would suffer the same fate if we took their infall 
motion towards larger objects into account. This is consistent 
with our observation that high-mass halos (which subsume the 
lower mass halos in their outskirts) undergo some physical ac- 
cretion in addition to the pseudo-evolution. 

A naive comparison of the mass function differences be- 
tween the forward and the backward evolution case at the low- 
mass end (see the bottom panels of Fig. [7]) seems to suggest 
that the effect of removing subhalos is very small. However, 
we note that for the case of forward evolution, it is extremely 
important to remove such subhalos. Due to their proximity to 
a larger host halo, the outskirts of their density profiles con- 
tain significant mass contributions from the host halo. As the 
virial radius increases towards lower redshift, such subhalos 
gain a significant fraction of the host halo mass. If they are 
not removed, this unphysical mass evolution results in a large 
scatter in the mass evolution histories of the low mass ha- 
los. This scatter can have a large eff ect o n the estimated mass 
function due to Eddingtonbias (see {33.21 i. Thus, not removing 
the 14% of halos which become subhalos by z — can lead 
to residuals of over 100% when comparing to the T08 mass 
function for the case of forward evolution. 

4. DISCUSSION 

In the past two sections, we have demonstrated that pseudo- 
evolution due to changing reference density has a significant 
impact on the overall evolution of SO mass (often called mass 
accretion history). In this section, we expand on some of the 
implications of this result for our understanding of the scaling 
relations between various observables and halo mass. This 
includes the concentration-mass relation, the relation between 
stellar content and halo mass, and scaling relations for galaxy 
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Fig. 8. — Evolution of the concentration-mass relation with redshift expected 
if the physical density around density peaks is unchanged over time. The solid 
lines show the analytical prediction of Equations {6}-0. From top to bottom, 
the lines correspond to redshifts to 1 in steps of 0.2. The dashed lines show 
the concentration-mass relation as a function of redshift from the physical 
model of Z09. This figure demonstrates that at low masses a large fraction of 
the observed evolution in the c—M relation is simply due to pseudo-evolution. 

clusters. 

4. 1 . The concentration-mass relation 

In the presence of pseudo-evolution, the virial radius of a 
halo grows with time, even though the halo's physical den- 
sity profile (and thus its scale radius, r s ) remain unchanged. 
However, since the virial radius does change due to pseudo- 
evolution, the concentration, c, grows at the same rate as the 
virial radius. Thus, we expect a significant evolution in the 
concentration-mass relation (hereafter c — M relation), even if 
the physical density profile of a halo remains unchanged. 

We use the static halo model to estimate the magnitude of 
this evolution. The prediction of the static halo model for two 
mass definitions is shown with solid lines Fig. [8] and is in 
qualitative agreement with the evolution observed in numeri- 
cal simulations such that the concentration of halos of a given 
mass decreases with increasing redshift. In reality, we expect 
that halos undergo some true physical evolution of mass due 
to physical accretion and merging, especially at the high-mass 
end, which will result in quantitative discrepancies between 
the true evolution of the c — M relation observed in numerical 
simulations and our static halo model predictions. 

The dashed lines in Fig.|8]show the redshift-dependent c—M 
relation obtained from the models of Z09 which have been 
calibrated to reproduce the evolution of this relation in nu- 
merical simulations. The comparison clearly shows that most 
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of the evolution in the c - M relation at the low-mass end can 
be accounted for by pseudo-evolution of halo radius at differ- 
ent redshifts. For high-mass halos, however, the Z09 model 
captures the c - M evolution due to their significant physical 
mass accretion, and is thus not completely reproduced by our 
model of pseudo-evolution. Note that the static halo model 
quantitatively reproduces the e volution of the c — M relation 
for low-mass halos obtained bv lBullock et al.l (120011) . 

4.2. Implications for galaxy formation 

The total stellar content (or the stellar light) that we observe 
as galaxies in a halo is the integral result of the complex inter- 
play between a variety of processes such as star formation, 
feedback from young stars, supernovae and supermassive 
black holes, and galactic outflows, all of which occur within 
dark matter halos. The redshift evolution of the scaling rela- 
tion between this stellar mass (or luminosity) and halo mass 
can provide important observational clues regarding these dif- 
ferent physical processes, in particular their efficiency as a 
function of halo mass. A number of studies have investigated 
the scaling relations between stellar mass and the mass of the 
halo t hey inhabit, and how these scaling relations change with 
time ([Conrov et alJ 120071: iBrown et alJl2008t IBehroozi et alJ 
| 2010tlMoster et al.l2010tlAbbas et alj2010t lWake et al.1201 li 
ILeauthaud et al.l 120121 : lYang et al.l l2012bl) . However, when 
connecting the observed evolution of the scaling relations to 
the underlying physics, it is crucial to account for the pseudo- 
evolution of halo mass. 

One of the striking implications of our results is that almost 
all of the mass evolution of galactic-sized halos (M2oop[z = 
0] < 10 12 /t'Mq) since z = 1 can be attributed to pseudo- 
evolution (see Fig. [3). The density profiles around the peaks 
of such halos have stayed static and not evolved physically 
(see a lso iPrada et aT1l2006h iDiemand et al.ll2007t ICuesta et alJ 
l2008h . Since physical accretion plays a minor role compared 
to pseudo-evolution, the impact of pseudo-evolution must be 
considered while interpreting the evolution of scaling rela- 
tions and relating them to the underlying physical processes. 

For example, the ratio of stellar mass to halo mass (SHMR), 
and its evolution with redshift, give a quantitative measure 
of how the star formation efficiency in a halo of given mass 
evolves with redshift. The peak of the star formation ef- 
ficiency lies at roughly Mawpfe = 0] - 10 12 /i~'M , and 
ha s been observed to shift to higher values from z — to 
1 dMoster et al.ll2010t IBehroozi et alll2010t ILeauthaud et all 
I2012D . However, equal halo masses at two different redshifts 
correspond to different physical density peaks due to pseudo- 
evolution. Thus, the rate at which the similar density peaks 
become inefficient can differ from the estimates at fixed halo 
mass. 

As a second example, let us consider the evolution of the 
SHMR at the more massive end. The stellar mass in such ha- 
los is dominated by satellite galaxies. Pseudo-evolution will 
lead to a constant fractional increase in both the stellar content 
and the halo mass if the distribution of satellite galaxies in and 
around halos, to first order, follows the matter density distri- 
bution, and if there is no radial segregation in the stellar mass 
of satellite galaxies. In this case, stellar and halo mass grow 
by the same factor, and the SHMR at the massive end un- 
dergoes an evolution with redshift w hich is qualitatively ver y 
similar to the evolution observed by ILeauthaud et ail ( 120121) . 
We will perform a quantitative comparison between the evo- 
lution of the SHMR due to pseudo-evolution of halo mass and 



the observed SHMR evolution in future work. 

4.3. Implications for cluster scaling relations 

As discussed in §Q] scaling relations between the baryonic 
properties of clusters, such as X-ray temperature, gas mass, or 
entropy, and the mass of the cluster's dark matter halo, are key 
to our understanding of clusters and their use in cosmology. 
The simplest model for these scaling relations relies on the 
assumption th at cluster halos collapse in a self-similar fashion 
dKaiserlll986l) . In this model, the temperature T, for example, 
scales with halo mass as 



where M denotes the mass within the radius R, and T is mea- 
sured at R (KB 12) . If the halo mass is defined as a spherical 
overdensity mass (mass definitions such as Msoop c or Afeoop, 
are commonly used for clusters), the above scaling relation 
can be expressed as 

T <* (A cPc ) 1/3 M 2/3 . (13) 

Noting that A c is a constant, but that the critical density 
evolves with E 2 (z), the evolution with redshift can be incor- 
porated into the scaling relation as 

T cc [E(z)M A ] 2 ' 3 . (14) 

Unfortunately, the E 2 (z) factor only accounts for the evolution 
of R (due to the evolution of the p c factor in Equation (fT3l), 
but not the pseudo-evolution of Ma. Thus, for a halo whose 
density profile remains constant, the scaling relation predicts 
that the temperature will increase with time without any par- 
ticular physical reason. 

As shown in Fig. [4] the mass evolution of Msoo Pc since 
Z = 1 is only partly due to pseudo-evolution, but also con- 
tains a large contribution from actual accretion. Nevertheless, 
pseudo-evolution accounts for a factor of almost two. If this 
contribution is not taken into account, pseudo-evolution will 
masquerade as a deviation from self-similar behavior. We in- 
tend to further investigate the impact pseudo-evolution has on 
scaling relations in the future. However, finding a mass def- 
inition which allows us to disentangle the effects of pseudo- 
evolution on cluster scaling relations is certainly a challeng- 
ing task, especially because different cluster observables have 
different dependencies on the exact boundary used to define a 
halo. 

5. CONCLUSION 

Several authors have pointed out that spherical over- 
density masses un dergo an evolution due to the changing 
reference density dPiemand et al.l 120071 : ICuesta et alJ 120081 : 
iRravtsov & Borganill2012l) " We have studied this spurious 
mass evolution quantitatively, and called it pseudo-evolution. 
Our main results are as follows. 

1 . We have demonstrated that a significant fraction of the 
halo mass growth since z = 1 is due to pseudo-evolution 
rather than actual physical accretion of matter. For low- 
mass halos ( < 10 12 /t'Mq), the pseudo-evolution ac- 
counts for almost all of the evolution in mass since 
z= 1. 

2. We found consistency between results derived from the 
simulated density profiles at z — and z = 1, indicating 
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that the physical shape of those profiles undergoes very 
little change over this epoch. We demonstrated that, at 
Z = 1, the infall velocities around low-mass halos are 
insufficient to facilitate significant growth, which di- 
rectly supports the observation of non-evolution of halo 



3. We have shown that a simple analytical model repro- 
duces the pseudo-evolution observed in the Bolshoi 
simulation to a few percent accuracy. 

4. We investigated the effect of pseudo-evolution on the 
halo mass function dN/dln(M), and found it to sim- 
ply shift the function toward higher masses. We pro- 
pose that the non-evolution of the mass function at low 
masses since z = 1 constitutes a fortuitous cancellation 
between pseudo-evolution and the absorption of small 
halos into larger halos. 

We have left some promising avenues of research for future 
investigations. For example, in this paper we restricted our- 
selves to halos with spherical overdensity definitions. How- 
ever, another popular way to identify and define halos is to 
use the friends-of-friends (FoF) algorithm, which relies on a 
linking length (which is fixed relative to the average inter- 
particle comoving separation) rather than an overdensity to 
define masses. It is well known that the density of FOF ha- 
los at their boundary depends upon linking length (s ee e.g. 
iFrenk et al.l 19881 iLukic et al.l2009tlMore et al.l201 lal) . Since 
the linking length parameter is constant in comoving coordi- 



nates, its physical length increases with time as the scale fac- 
tor (a = [1 + z] _1 ). This implies that for a static halo density 
profile, the extent of the FoF halo will increase with time, 
leading to pseudo-evolution . While explaining the mass evo- 
lution history of FOF halos. iFakhouri & Mai (1201 Oh disentan- 
gle the growth into accretion of resolved halos and the accre- 
tion of a diffuse component. It is clear that the diffuse compo- 
nent will include a contribution from pseudo-evolution, thus 
overestimating the fraction of actually accreted diffuse matter 
compared to the fraction accreted from merging halos. There- 
fore, the pseudo-evolution of FoF halos will need to be care- 
fully investigated. 

The eventual goal of such investigations will be to find a 
mass definition which does not suffer from pseudo-evolution, 
which can be measured observationally and in simulations, 
and which allows the formulation of meaningful scaling rela- 
tions. We hope to address this question in future work. 
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APPENDIX 
MASS EVOLUTION 

The halo mass at a given redshift zo is given by the integration of the density profile from the center to the radius of the halo 
R&, which depends upon the definition of the halo, 

/~«a(zo) 

M A (z ) = p(r',z)4nr'W . (Al) 

Jo 

The halo mass evolution from redshift zo to a given redshift z can be decomposed into the following components 

/~«'a(z) /~«'a(z) /-fiAfe) 

O) = M (zo) + p(r ,zo)4n r ' 2 dr' + [p(r' ,z) - p(r' ,z )]4nr' 2 dr' + p(r' ,z)4nr' 2 dr' (A2) 

J«a(zo) Jo Jr'a(z) 

where R' A represents the estimate of the radius at redshift z calculated based upon the static halo model, which is slightly different 
from the true radius at redshift z due to physical evolution of the halo. The first integral represents the evolution of mass just 
due to the changing mass definition. The second integral accounts for the change in the halo mass within the radius R'a(z) due 
to change in density from redshift zq to z- The third integral just adds up the mass between the halo boundary predicted by the 
static halo model, R' A and the true boundary, R&. If halos are entirely static as assumed in our analytical model, only the first 
integral contributes to the mass evolution and the rest of the integral terms are exactly zero. This prediction from the static halo 
model is shown as solid lines in Fig. [3] while the total mass evolution from models of Z09 is shown using blue dashed lines. The 
comparison shows that at low halo masses, the mass evolution is indeed dominated by the first integral term. 



